Observer dependence for the phonon content of the sound field hving on the eff'ective 
curved space-time background of a Bose-Einstein condensate 
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We demonstrate that the ambiguity of the particle content for quantum fields in a generally curved 
space-time can be experimentally investigated in an ultracold gas of atoms forming a Bose-Einstein 
condensate. We explicitly evaluate the response of a suitable condensed matter detector, an "Atomic 
Quantum Dot," which can be tuned to measure time intervals associated to different effective acous- 
tic space-times. It is found that the detector response related to laboratory, "adiabatic," and de 
Sitter time intervals is finite in time and nonstationary, vanishing, and thermal, respectively. 
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I. INTRODUCTION 

The feat of reproducing certain features of the physics 
of relativistic quantum fields on curved space-time back- 
grounds P, 01 , in the laboratory, now seems closer than 
ever before. Due to the realization of the fact that 
phonons or generally "relativistic" quasiparticles, prop- 
agating in a nonrelativistic background fluid, experience 
an effective curved space-time [3,SlEIEI3' various exotic 
phenomena of the physics of classical as well as quantum 
fields in the curved space-times of gravity are getting 
within the reach of laboratory scale experiments • 

The particle content of a quantum field state in curved 
or flat space-time depends on the motional state of the 
observer. A manifestation of this observer dependence is 
the Unruh-Davies effect, which consists in the fact that a 
constantly accelerated detector in the Minkowski vacuum 
responds as if it were placed in a thermal bath with tem- 
perature proportional to its acceleration This effect 
has eluded observation so far: The value of the Unruh 
temperature Tunruh = [ft/(27rA:BCL)]a = 4K x a[lO^°50] , 
where a is the acceleration of the detector in Minkowski 
space (50 is the gravity acceleration on the surface of the 
Earth), and cl the speed of light, makes it obvious that 
an observation of the effect is a less than trivial under- 
taking. Proposals for a measurement with ultraintense 
short pulses of electromagnetic radiation have been put 
forward in, e.g., Refs. [9l[l0|. 

In the following, we shall argue that the observer 
dependence of the particle content of a quantum field 
state in curved space-time, related to the familiar non- 
uniqueness of canonical field quantization in Riemannian 
spaces can be experimentally demonstrated in the 
readily available ultralow temperature condensed matter 
system Bose-Einstein condensate (BEC). More specifi- 
cally, we argue that the Gibbons-Hawking effect a 
curved space-time generalization of the Unruh-Davies ef- 
fect, can be observed in an expanding BEC. To explicitly 
show that particle detection depends on the detector set- 
ting, we construct a condensed matter detector tuned to 



time intervals in effective laboratory and de Sitter space- 
times, respectively 0]. We show that the detector re- 
sponse is strongly different in the two cases, and associ- 
ated to the corresponding effective space-time. Further- 
more, we describe a system of space-time coordinates in 
which there is no particle detection whatsoever taking 
place, which we will call the "adiabatic" basis, and which 
simply corresponds to a detector at rest in the (confor- 
mal) Minkowski vacuum. 



II. HYDRODYNAMICS IN AN EXPANDING 
CIGAR-SHAPED BEC 

It recently became apparent that among the most suit- 
able systems for the simulation of quantum phenomena 
in effective space-times are Bose-Einstein condensates 
(BECs) El 111 EE m m. They offer the primary 
advantages of dissipation-free superflow, high controlla- 
bility, with atomic precision, of the physical parameters 
involved, and the accessibility of ultralow temperatures 
|l9l l20l | . Even more importantly, in the context under 
consideration here, the theory of phononic quasiparticles 
in the inhomogeneous BEC is kinematically identical to 
that of a massless scalar field propagating on the back- 
ground of curved space-time in D + 1 dimensions. 

It was shown by Unruh that the action of the phase 
fluctuations $ in a moving inhomogeneous superfluid 
may be written in the form d (we set h 
where m is the mass of a superfluid constituent particle): 
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(1) 



Here, v{x,t) is the superfluid background velocity, 
c{x,t) — y^gpo{x, t) is the velocity of sound, where g 
is a constant describing the interaction between the con- 
stituent particles in the superfluid {1/g is the compress- 
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ibility of the fluid), and pQ{x,t) is the background den- 
sity. In the second line of ^ , the conventional hydrody- 
namic action is identified with the action of a minimally 
coupled scalar field in an effective curved space-time. 
Furthermore, the velocity potential of the sound pertur- 
bations in the BEC satisfies the canonical commutation 
relations of a relativistic scalar field We therefore 

have the exact mapping, on the level of kinematics, of the 
equation of motion for phononic quasiparticles in a non- 
relativistic superfluid, to quantized massless scalar fields 
propagating on a curved space-time background with lo- 
cal Lorentz invariance. 

To model curved space-times, we will consider the evo- 
lution of the BEC if we change the harmonic trapping 
frequencies with time. For a description of the expansion 
(or contraction) of a BEC, the so-called scaling solution 
approach is conventionally used One starts from 

a cigar-shaped BEC containing a large number of con- 
stituent particles, i.e., which is in the so-called Thomas- 
Fermi (TF) hmit H^. According to [111, the evolution 
of a Bose-condensed atom cloud under temporal varia- 
tion of the trapping frequencies uj\\(t) and uj±{t) (in the 
axial and radial directions, respectively) can then be de- 
scribed by the following solution for the condensate wave 
function 



^ = 



TF 



b^Vb 



exp 



bz"^ b^T^ 
gpoix = 0, t)dt + i— + 



(2) 

Here, b± and b are the scaling parameters describing the 
condensate evolution in the radial (f) and axial (z) di- 
rections, cf. Fig.^ The initial {b = b± = 1) mean-field 
condensate density is given by the usual TF expression 



|*tf|^ = Ptf('', z) = pm{l 



(3) 



Here, pm is the maximum density (in the center of the 
cloud) and the squared initial TF radii are i?^ = 2/i/w| 

and R\ = The initial chemical potential /i = 

Pm5, where g — 47ras, and as is the scattering length 
characterizing atomic collisions in the (dilute) BEC. In 
our cylindrical 3D trap, we have for the initial central 
density 
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The condition that the TF condition be valid implies 
that /i ^ wii,^^. The solution ^ of the Gross- 
Pitaevskiimean-field equations becomes exact in this TF 
limit, independent of the ratio uj\\Ilo\^. However, the so- 
lution becomes exact also in the limit that 0, 
independent of the validity of the TF limit, the system 
then acquiring an effectively two-dimensional character 
[25I l2^ . We will see below that in the opposite limit 
of a;||/w_L there is an "adiabatic basis" in which no 
axial excitations are created during the expansion, i.e.. 
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FIG. 1: Expansion of a cigar-shaped Bose-Einstein conden- 
sate. The stationary horizon surfaces are located at ±zh, 
respectively. The thick dark lines represent lasers creating an 
optical potential well in the center. 



with respect to that basis there are, in particular, no 
unstable solutions possible, implying the stability of the 
expanding gas against perturbations. 

According to lO, the condensate density evolves as 



PQ{r,z,t) = 



PTF {rVbl,zyb^) 



and the superfluid velocity 



V = rEr 

b I 
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(4) 



(5) 



is the gradient of the condensate phase in Eq. (O. It 
increases linearly to the axial and radial boundaries of 
the condensate. 

The excitations in the limit aju/cij^ — s- were first stud- 
ied in [Slj. The description of the modes is based on an 
adiabatic separation for the axial and longitudinal vari- 
ables of the phase fluctuation field: 



^{r,z,t) = ^(?!)„(r)x„(2;,t), 



(6) 



where (pnii") is the radial wavcfunction characterized by 
the quantum number n (we consider only zero angular 
momentum modes). The above ansatz incorporates the 
fact that for strongly elongated traps the dynamics of the 
condensate motion separates into a fast radial motion and 
a slow axial motion, which are essentially independent. 
The Xn(z, i) are the mode functions for travelling wave 
solutions in the z direction (plane waves for a condensate 
at rest read Xn exp[— ie„ /ji -t- kz]). The radial motion 
is assumed to be "stiff" such that the radial part is ef- 
fectively time independent, because the radial time scale 
for adjustment of the density distribution after a pertur- 
bation is much less than the axial oscillation time scales 
of interest. The ansatz © works independent from the 
ratio of healing length and radial size of the BEC cigar. 
In the limit that the healing length is much less than the 
radial size, TF wave functions are used, in the opposite 
limit, a Gaussian ansatz for the radial part of the wave 
function 0„(r) is appropriate. 
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For axial excitations characterized by a wavelength A = 
27r/fc exceeding the radial size R± of the condensate, we 
have kR± <C 1, and the dispersion relation reads, in the 
TF limit for the radial wave function 123 



- 2cjin{n + l) + ^{kR^y 
^ 2Luln{n + 1) + clk^ , 



(7) 



where cq — \/ n/2. Observe that the central speed 
of sound Co is reduced by a factor \/2 from the well- 
established value y/gpm ~ \/m for infinitely extended 
liquid 

The Eqs. JHl and Q can be generalized for an ex- 
panding condensate. Substituting in Eq. ^ the rescaled 
radial wavcfunction 0„ = 0n(''/&±), and inserting the 
result into the action integrating over the radial co- 
ordinates, we find the following effective action for the 
axial modes of a given radial quantum number n: 



dtdz 



23 



+Ml{z)xl] , (8) 



where the common "conformal" factor Cn{z) is given by 



J r<rm 



(9) 



The integration limits are fixed by the z dependent radial 
size of the cigar = ~ The averaged 

speed of sound reads 



cl{z,t) 
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r /)2 

^7i^j_ Jr<rr 



d rpo(j)„ 



(10) 



and the (space and time dependent) effective mass term 
is, for a given radial mode, obtained to be 



III. THE 1+lD DE SITTER METRIC IN THE 
CONDENSATE CENTER 

We identify the action with the action of a mini- 
mally coupled scalar field in 1-l-lD, according to Eq. 
Remarkably, such an identification is possible only close 
to the center z — 0. The contravariant 1+lD metric may 
generally be written as Q 



-1 



(14) 



where Ac is some arbitrary (space and time dependent) 
factor and c = c(z = 0). Inverting this expression to get 
the covariant metric, we obtain 



^fiiy — Aq 



-V, 1 



(15) 



The term -v/^g^" contained in the action ^ gives the 
familiar conformal invariance in a 1-l-lD space-time, i.e., 
the conformal factor Ac drops out from the action and 
thus does not influence the classical equations of motion. 
We therefore leave out Ac in the formulae to follow, but 
it needs to be borne in mind that the metric elements are 
defined always up to the factor Ac- 

The actions and ^ can be made consistent if we 
renormalize the phase field according to $ = Z$ and 
require that 
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1 



(16) 



holds. The factor Z docs not influence the equation of 
motion, but does influence the response of a detector (see 
section HV Al belowl . In other words, it renormalizes the 
coupling of our "relativistic held" $ to the lab frame 
detector. More explicitly, Eq. \i<6\ leads to 



d^rpQ [dr(f>n 



(11) 



The phonon branch of the excitations corresponds to the 
gapless n = solution of Eq. |7J . In this case the radial 
wavefunction (f>Q does not depend on the radial variable 
r [231, and the mass term vanishes, Mq = 0. We then 
obtain the following expressions. 



Co(z) = 1 



b^Rj 



and for the z dependent speed of sound (c = Cn=o)'- 
bib \ b^Rj J ■ 



c^{z,t) 



(12) 



(13) 



We will see below that we need these expressions in the 
limit z only, because only in this limit we obtain the 
exact mapping of the phonon field to a quantum field 
propagating in a 1-t-lD curved space-time. 



Vb 
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2Z^^ 



UJJ_ 



— = B ^ const. 



(17) 



According to the above relation, we have to impose that 
the expansion of the cigar in the perpendicular direction 
proceeds like the square root of the expansion in the axial 
direction. The constant quantity B can be fixed exter- 
nally (by the experimentalist), choosing the expansion of 
the cloud appropriately by adjusting the time dependence 
of the trapping frequencies W|| (t) and U!±(t), according to 
the scaling equations 0i 1^ 



b + Luj{t)b = 



5253' 

,2 



LO 



5355/2 ■ 



(18) 



Since both Cq and c depend on z, an effective space- 
time metric for the axial phonons can be obtained only 
close to the center of the cigar-shaped condensate cloud. 
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This is related to our averaging over the physical per- 
pendicular direction, and does not arise if the excitations 
are considered in the full D-dimensional situation, where 
this identification is possible globally. Also note that the 
action (jSJ does not contain a curvature scalar contribu- 
tion of the form oc XnV^& that it only 
possesses trivial conformal invariance |29j |. 

We now impose, in addition, the requirement that the 
metric is identical to that of a 1-l-lD universe, with a met- 
ric of the form of the de Sitter metric in 3-l-lD [T^ |2]|. 
We first apply the transformation codt — c{t)dt to the line 
element defined by (|15|l . connecting the laboratory time t 
to the time variable i. Defining Vz/c = ^/Az = {Bb/co)z 
(note that the "dot" on b and other quantities always 
refers to ordinary laboratory time), this results, up to 
the conformal factor A^, in the line element 



ds^ 



-cl{l - Az^)dP - 2cozVXdidz + dz'^. (19) 

We then apply a second transformation codr — c^dt + 
z^/Adz/{l — Az^), with a constant A. We are thus led to 
the 1+lD de Sitter metric in the form |0 



ds2 ^ -cl (1 - Az2) dr' + (1 - Az2) 



dz^ . (20) 

The transformation between t and the de Sitter time r 
(on a constant z detector, such that dt = dt), is given by 



— = exp[i35T] 



(21) 



where the unit of lab time tt^ ^ Wn ^ is set by the initial 
conditions for the scaling variables b and b±. The tem- 
perature associated with the effective metric (I20II is the 
Gibbons-Hawking temperature 



dS 



2tt 



(22) 



The "surface gravity" on the horizon has the value an = 
CqVA = CoBb, and the stationary horizon(s) are located 
at the constant values of the z coordinate 



±Zh = ±i?|| 



2^A 



(23) 



Combining (|22() and (|23|l . we see that zb_/R\\ is small 
if wii/Tds ^ 47r. Therefore, the de Sitter temperature 
needs to be at least of the order of wn for the horizon lo- 
cation(s) to be well inside the cloud. The latter condition 
then justifies neglecting the z dependence in Co and c in 
Eq. Hl()|l . Though there is no metric "behind" the hori- 
zon, i.e. at large z, this should not affect the low-energy 
behavior of the quantum vacuum "outside" the horizon. 



IV. DETERMINING THE PARTICLE 
CONTENT OF THE QUANTUM FIELD 

The particle content of a quantum field state depends 
on the observer P, 0, 0. To detect the Gibbons- 
Hawking effect in de Sitter space, one has to set up a 



detector which measures frequencies in units of the in- 
verse de Sitter time r, rather than in units of the in- 
verse laboratory time t. We will show that the de Sit- 
ter time interval dr = dt/bB ~ dt/\/bb± can be effec- 
tively measured by an "Atomic Quantum Dot" (AQD) 
|22ll23i] . The measured quanta can then, and only then, 
be accurately interpreted to be particles coming from a 
Gibbons-Hawking type process with a constant de Sitter 
temperature (|22|l . We then use the tunability to other 
time intervals feasible with our detector scheme, and con- 
trast this de Sitter result with what the detector "sees" 
if tuned to laboratory and "adiabatic" time intervals. 

The AQD can be made in a gas of atoms possessing 
two hyperfine ground states a and /3. The atoms in the 
state a represent the superfluid cigar, and are used to 
model the expanding universe. The AQD itself is formed 
by trapping atoms in the state /3 in a tightly confining op- 
tical potential Vopt- The interaction of atoms in the two 
internal levels is described by a set of coupling parame- 
ters gcd = '^T^O'cd (c, d — {a, /?}), where acd are the s-wave 
scattering lengths characterizing short-range intra- and 
inter-species collisions; Qaa = g, CLaa = a^. The on-site 
repulsion between the atoms (3 in the dot is U gpp/l^, 
where / is the characteristic size of the ground state wave- 
function of atoms P localized in Vopt- In the following, 
we consider the collisional blockade limit of large J7 > 0, 
where only one atom of type f3 can be trapped in the dot. 
This assumes that U is larger than all other relevant fre- 
quency scales in the dynamics of both the AQD and the 
expanding superfluid. As a result, the collective coordi- 
nate of the AQD is modeled by a pseudo-spin- 1/2, with 
spin-up/spin-down state corresponding to occupation by 
a single/no atom in state (3. 

We first describe the AQD response to the condensate 
fiuctuations in the Lagrangian formalism, most familiar 
in a field theoretical context. The detector Lagrangian 
takes the form 



^AQD = i{—'n*]rj-[-A + ga[3{po{0,t) + Sp)]ri*r] 



exp ~i / gpo{0,t')dt' + ii 



■ exp 



i I gpoiO,t')dt' - iS(b 71 ]. (24) 



Here, A is the detuning of the laser light from resonance, 
Pa(z = 0,i) is the central mean- field part of the bath 
density, and I is the size of the AQD ground state wave 
function. The detector variable 77 is an anticommuting 
Grassmann variable representing the effective spin de- 
gree of freedom of the AQD. The second and third lines 
represent the coupling of the AQD to the surrounding 
superfluid, where Scf) and 6p are the fluctuating parts of 
the condensate phase and density at z = 0, respectively. 
The laser intensity and the effective transition matrix el- 
ement combine into the Rabi frequency below we will 
make use of the fact that can easily experimentally be 
changed as a function of laboratory time t, by changing 
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the laser intensity with t. 

To simplify H24|l . we use the canonical transformation 



V 



T] exp 



— i / gpo{0,t')dt' + iS(j) 





(25) 



The above transformation amounts to absorbing the su- 
perfluid's chemical potential and the fluctuating phase S(f> 
into the wave function of the AQD, and does not change 
the occupation numbers of the two AQD states. The 
transformation (|25|l gives the detector Lagrangian the 
form 



^AQD 



^ ^ -fj* ) fj - ny/Mo^ iv + vl (26) 



-A + {gai3 - g)pa{0, t) + gapSp + 



dt 



77 77. 



The laser coupling (second term in the first line) scales 
as and hence like the de Sitter time interval in 

units of the laboratory time interval, dr/dt. We suggest 
to operate the detector at the time dependent detuning 
A(i) = (ffa/3 - <7)Po(0,t) = [g^p - g)p^l[h'^B^e), which 
then leads to a vanishing of the first two terms in the 
square brackets of (j^ . 

We now introduce back the wave function of the AQD 
stemming from a Hamiltonian formulation, ip = "0/31/3) + 
'ipa\oi). An "effective Rabi frequency" may be defined to 
be uiQ — 2Q,yJ Pml^', at the detuning compensated point, 
we then obtain a simple set of coupled equations for the 
AQD amplitudes 

= Y^" + ^^^'^' '^^^^''^ (^^^ 

where r is the de Sitter time. 

We have thus shown that the detector Eqs. H27|l are 
natural evolution equations in de Sitter time r, if the 
Rabi frequency is chosen to be a constant, indepen- 
dent of laboratory time t. We will see in section IIV Bl 
that, adjusting in a certain time dependent manner, 
within the same detector scheme, we can reproduce time 
intervals associated to various other effective space-times. 

The coupling of the AQD to fluctuations in the super- 
fluid is described by the potential 



5V{T)^{g^p-g)Bb{T)5p{T). 



(28) 



Neglecting the fluctuations in the superfluid, the level 
separation implied by H27() is f^nd the eigenfunctions 
of the dressed two level system are |±) = {\a) ± |/3))/-\/2. 
The quantity loq therefore plays the role of a frequency 
standard of the detector. By adjusting the value of the 
laser intensity, one can change ojq, and therefore probe 
the response of the detector for various phonon frequen- 
cies. Note that if gap is very close to g, to obtain the 
correct perturbation potential, higher order terms in the 
density fluctuations have to be taken into account in the 
Rabi term of if^ . 



To describe the detector response, we first have to 
solve the equations of motion ||SJ) for the phase fluctu- 
ations, and then evaluate the conjugate density fluctua- 
tions. The equation of motion 5Sq/5xo = is, for time 
independent S, given by 



(29) 

where — z/b is the scaling coordinate. Apart from the 
factor 6*0(2:6), stemming from averaging over the perpen- 
dicular direction, this equation corresponds to the hy- 
drodynamic equation of phase fluctuations in inhomoge- 
neous superfluids |30j |. At i — > — cx), the condensate is in 
equilibrium and the quantum vacuum phase fluctuations 
close to the center of the condensate can be written in 
the following form 



Xo = 



9 



4Co(0)i?||eo,fc 



dk exp [— ico.fci + ikz] + H.c, (30) 



where afc,aj, are the annihilation and creation operators 
of a phonon. The intial quantum state of phonons is the 
ground state of the superfluid and is annihilated by the 
operators dk- With these initial conditions, the solution 
of (El is 



Xo 



4Co(0)i?||eo,fc 
-l-H.c. 



dk exp 



ikzh 



(31) 



The solution for the canonically conjugate density fluc- 
tuations, consequently, is 



5p 



d 



ACo(,0)Riigdt 
-hH.c. 



ftfe exp 



dt'ep^f. 
Bb^ 



ikzh 



(32) 



The Eqs. H31|l and H32() completely characterize the n = 
evolution of the condensate fluctuations. Observe that 
the evolution proceeds without frequency mixing in the 
adiabatic time interval deflned by dr^, = dtjBb'^ (the 
"scaling time" interval dt/B^b^ deflned in jsj is pro- 
portional to this adiabatic time interval). Therefore, in 
the "adiabatic basis," no frequency mixing occurs and 
thus no quasiparticle excitations are created. This hints 
at a hidden (low energy) symmetry, in analogy to the 
(exact) 2+lD Lorentz group S0(2,l) for an isotropically 
expanding BEG disk, discussed in |25| . 



A. Detection in de Sitter time 

The coupling operator SV causes transitions between 
the dressed detector states |-|-) and |— ) and thus can 
be used to effectively measure the quantum state of the 
phonons. We consider the detector response to fluctu- 
ations of by going beyond mean-fleld and using a 
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perturbation theory in SV. There are two physically 
different situations. The detector is either at t = in 
its ground state, (|a) + |/3))/\/2, or in its excited state, 
(|a) - \f3))/V2. We define P+ and P_ to be the prob- 
abilities that at late times t the detector is excited re- 
spectively de-excited. Using second order perturbation 
theory in SV, we find that the transition probabilities for 
the detector may be written 

gepM I gap 1 \ d2 

^4i?j|Co(0) V 5 

where the absolute square of the transition matrix ele- 
ment is given by 



P± 



— -1 B'\T±\\ (33) 



dr 



■ exp 



dr' 

W) 



Calculating the integrals, we obtain 



P± = J 



9a() 

9 



2BhR\\Co{Q) 
where the (formally divergent) sum 



ILUqT 



71b 

1 + n-B 



and the factors 



1 



(34) 



,(35) 



(36) 



(37) 



exp[a;o/Tds] 

are Bose distribution functions at the de Sitter temper- 
ature H22|) . We conclude that an expansion of the con- 
densate in z direction, with a constant rate faster than 
the harmonic trap oscillation frequency in that direction, 
gives an effective space-time characterized by the de Sit- 
ter temperature Tds- 

We now show that J is proportional to the total de 
Sitter time of observation, so that the probability per 
unit time is a finite quantity [s^ . At late times, the 
detector measures phonon quanta coming, relative to its 
space-time perspective, from close to the horizon, at a 
distance Sz = zji — z <$:. zn = A~^/^. The trajectory of 
such a phonon in the coordinates of the de Sitter metric 
(|20|l . at late times r, is given by (cf. Fig.|2J) 

In [^1 = 2VAcot. (38) 
L dzi 

This implies that the central AQD detector measures 
quanta which originated at the horizon with the large 
shifted frequency 



eo,fc — 



V2Al/4^;,l/2 ^ 



exp[coTV A 



(39) 



Making use of the above equation, we rewrite the sum- 
mation over k in (|36|l as an integral over detector time: 



^0 



■KCO 

R\\i^o 

TTCO 



TTCO 



-\/Aco / dr 



Bhr. 



(40) 




FIG. 2: Typical trajectories of phonons {ds^ = 0) in the 
de Sitter metric 1201 . The path taken by phonon I, which 
is at early times propagating near the horizon, is described 
in the text. The path taken by phonon II, which approaches 
the horizon surface — 2h at late times r, does not lead to an 
excitation of the de Sitter detector placed at z = 0. 



Therefore, the probabilities per unit detector time (de 
Sitter time) read, where upper/lower entries refer to 
respectively: 



dP± 
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1 B 



2 5^0 



2Co(0)co i 1 + nB 



rtB 



(41) 



They are finite quantities in the limit that t — s- oo. In 
laboratory time, the transition probabilities evolve ac- 
cording to 



P±{t) 



Pn m 



TlB 
1 + TIB 



(42) 



where, from relation Pq = Z^[{ga,f3/g-l)B]'^/2. We 
see that the detector response is, as it should be, propor- 
tional to Z^, the square of the renormalization factor of 
the phase fluctuation field. 

The absorption and emission coefficients dP±/dT sat- 
isfy Einsteinian relations. Therefore, the detector ap- 
proaches thermal equilibrium at a temperature Tds on a 
time scale proportional to Z^^lj^^. Our de Sitter AQD 
detector thus measures a stationary thermal spectrum, 
even though its condensed matter background, with labo- 
ratory time t, is in a highly nonstationary motional state. 
Since Z^ oc ■\/ Pmof (w_l//x)^, not-too-dilute condensates 
with llJj_ ^ ^ (i.e., close to the quasi-lD regime [23|) are 
most suitable to observe the Gibbons-Hawking effect. 

The verification of the fact that a thermal detector 
state has been established proceeds by the fact that the 
two hyperfine states a and (3 are spectroscopically differ- 
ent states of the same atom, easily detectable by mod- 
ern quantum optical technology. When the optical po- 
tential is switched on, the atoms are in the empty a 
state originally, which is an equal-weight superposition 
of 1+) = (|a) + m/V2 and |-) = (|a) - |/3))/V2. The 
thermalization due to the Gibbons-Hawking effect takes 
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place in the dressed state basis consisting of the two de- 
tector states, i.e. of the states |+) and |— ), on a time scale 
given by the quantities P± in Eq. 1)42(1 . For the labora- 
tory observer, the Gibbons-Hawking thermal state will 
thus appear to cause damping of the Rabi oscillations on 
the thermalization timescale, i.e., friction on the coher- 
ent oscillating motion between the two detector states 
occurs, due to the thermal phonon bath perceived by the 
detector. The occupation of the detector states can be 
measured directly using atomic interferometry: A tt/2- 
pulsc brings one of them into the filled (/?) and the other 
into the empty (a) state. To increase the signal to noise 
ratio, one could conceive of manufacturing a small array 
of AQDs in a sufficiently large cigar-shaped host super- 
fluid, and monitor the total population of (3 atoms in this 
array. 



B. Detection in laboratory time 

We contrast the above calculation with the response 
the AQD detector would see if tuned to laboratory 
tim e. This can be realized if we let ft cx t, such that 
fly/ po{0, t) = n^/p^/{Bbt) = const., in the Rabi term 
on the right-hand side of (|26|) . is time independent. The 
detector has, therefore, dt as its natural time interval in 
this setting. The Painleve-GuUstrand metric (|15() in pure 
laboratory frame variables, assuming B^r ^ "j^^/^^ like 
in the derivation of the de Sitter metric H2U|I . reads 



ds^ 



1 - Az') dt' 



2z 

T 



-dzdt + dz 



(43) 



The metric (|43f) is asymptotically, for large t, becoming 
that of Galilei invariant ordinary ID laboratory space, i.e. 
is just measuring length along the z direction, because 
the speed of sound in the ever more dilute gas decreases 
like 1/t and the "phonon ether" becomes increasingly less 
stiff. 

The transition probabilities for absorption re- 
spectively emission are now given by P± ~ 
(5eo,fc/4i?||Co(0))(ga;3/g-l)'|f±|2, where the ma- 
trix elements are, cf. Eq. (|34|l . 



|T+|2 



dt 



exp 



±«eo,fc 



dt' 



iujot 



(44) 

Substituting the adiabatic time interval rfra = dt/ [Bb^t') 
leads for large t to 



(45) 



where Toa = dt/{Bb^t^). The transformation to adi- 
abatic time maps t e [— oo, -|-oo] onto Ta G [— oo, tq^] and, 
by further substituting y — eo,fc(''a — Toa), we have 



dy exp 



VT 



Bb^ y 



(46) 



The integral is a linear combination of Bessel functions. 
To test its convergence properties, we are specifically in- 
terested in the large eofe limit. Performing a stationary 
phase approximation for large A — zLujQeok / Bb'^ , we have 
for positive A (absorption) that the integral above be- 
comes J{A) — (ttvC?)^/^ exp[— 2\/C4] and for negative A 
(emission) J (A) = (7ry^|A|)^/^. The final result then is 



P+ = 



Bp 



/o V eofc-DO 



too ^1 exp 



Bb^ 




where Epi fi is the ultraviolet cutoff in the integral for 
the emission probability P_ , the "Planck" scale of the 
superfluid. Because of the convergence of the absorption 
integral for P+, the total number of particles detected 
remains finite, and there are no particles detected by the 
effective laboratory frame detector at late times. This is 
in contrast to the de Sitter detector, which according to 
(|41|l still detects particles, in a stationary thermal state. 

There is a detector setting which corresponds to a de- 
tector at rest in the Minkowski vacuum. This setting 
is represented by the adiabatic basis, with time interval 
defined by dra = dt/Bb', realizable with the AQD by 
setting the Rabi frequency f2 oc 1/t. Then, no particles 
whatsoever are detected, i.e., no frequency mixing of the 
positive and negative frequency parts of H32|l does take 
place. The associated space-time interval 



2 J 2 



(48) 



is simply that of (conformally) flat Minkowski space in 
the spatial scaling coordinate Zb and adiabatic time co- 
ordinate T„ . 



V. SUMMARY AND CONCLUSIONS 

We summarize the effective space-times considered in 
this article, and the associated time intervals in table 
^ The major observation of the present investigation 
is, that the physical nature of the effective space-time 
considered in the condensed matter system reflects itself 
directly in the quasiparticle (phonon) content measured 
by a detector which has a natural time interval equal to 
the time interval of this particular effective space-time. 
We have demonstrated that the notion of observer de- 
pendence can be made experimentally manifest by an 
Atomic Quantum Dot placed at the center of a linearly 
expanding cigar-shaped Bose-Einstein condensate, which 
has a tunable effective time interval. There is thus opened 
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Effective Space Time Interval Phonons Detected 

Laboratory dt Yes (Nonthermal) 

de Sitter dr = dt/Bb Yes (Thermal) 

Adiabatic dr^ = dt/Bb'^ No 

TABLE L Various time intervals effectively measuring labo- 
ratory, de Sitter, and adiabatic time in a I+ID BEC, respec- 
tively, where 6 = b{t), and B is independent of lab time. The 
third entry specifies if the AQD detector, tuned to the given 
time interval, detects phonons. In the laboratory frame, the 
detector has a small, nonthermal response for a finite amount 
of (initial) laboratory time, cf. Eq. 14711 . 

up a possibility to confirm experimentally, in an effective 
curved space-time setting, that, indeed, "A particle de- 
tector will react to states which have positive frequency 
with respect to the detector's proper time, not with re- 
spect to any universal time 0." 

The Gibbons-Hawking effect in the BEC is intrinsically 
quantum: The signal contains a "dimcnsionless Planck 
constant," i.e. the gaseous (loop expansion) parameter 
\/ Pmcfig- This implies that a reasonable signal-to-noise 
ratio can be achieved only by using initially dense clouds 



with strong interparticle interactions. On the other hand, 
the phononic quasiparticles of the superfluid can be re- 
g arded a s non-interacting only in a first approximation in 
yTWif ^ 1- The effects of self- interaction between the 
phonons, induced by larger values of the gaseous param- 
eter, can lead to decoherence and the relaxation of the 
phonon subsystem. The same line of reasoning applies 
to the evolution of quantum fields in the expanding uni- 
verse. The interactions between quasiparticle excitations 
and their connection to decoherence processes in cosmo- 
logical models of quantum field propagation and particle 
production is, therefore, an important topic for future 
work. 



Acknowledgments 

We acknowledge helpful discussions with R. Parentani 
and R. Schiitzhold. P. O.F. has been supported by 
the Austrian Science Foundation FWF and the Rus- 
sian Foundation for Basic Research RFRR, and U. R. F. 
by the FWF. They both gratefully acknowledge support 
from the ESF Programme "Cosmology in the Labora- 
tory." 



[1] W.G. Unruh, Phys. Rev. D 14, 870 (1976). 

[2] N.D. Birrell and P. C. W. Davies, Quantum Fields in 

Curved Space (Cambridge University Press, Cambridge, 

England, 1984). 
[3] W.G. Unruh, Phys. Rev. Lett. 46, 1351 (1981). 
[4] M. Visser, Class. Quantum Grav. 15, 1767 (1998). 
[5] R. Schiitzhold and W.G. Unruh, Phys. Rev. D 66, 

044019 (2002). 

[6] U. R. Fischer and M. Visser, Ann. Phys. (N.Y.) 304, 22 
(2003)); Phys. Rev. Lett. 88, 110201 (2002). 

[7] P. Nurowski, E. Schucking, and A. Trautman, "Rela- 
tivistic gravitational fields with close Newtonian analogs", 
Chapter 23 in "On Einstein's Path", Essays in Honor of 
Engelbert Schucking, Editor A. Harvey (Springer, New 
York, 1999). 

[8] G. E. Volovik, The Universe in a Helium, Droplet (Oxford 
University Press, Oxford, 2003); for a broad overview on 
effective space-times in condensed matter systems. 
[9] E. Yablonovitch, Phys. Rev. Lett. 62, 1742 (1989). 
[10] P. Chen and T. Tajima, Phys. Rev. Lett. 83, 256 (1999). 
[11] S.A. Fulling, Phys. Rev. D 7, 2850 (1973). 
[12] G.W. Gibbons and S.W. Hawking, Phys. Rev. D 15, 
2738 (1977). 

[13] We will use the notions of particle and quasiparticle inter- 
changeably, as far as excitations of quantum fields obey- 
ing relativistic respectively pseudo-relativistic wave equa- 
tions are concerned. The particles (atoms) forming the 
BEC will be called "constituent particles." 

[14] W.G. Unruh and R. Schiitzhold, Phys. Rev. D 68, 
024008 (2003). 

[15] C. Barcelo, S. Liberati, and M. Visser, Class. Quantum 
Grav. 18, 1137 (2001); Phys. Rev. A 68, 053613 (2003). 



[16] L.J. Garay, J.R. Anglin, J.I. Cirac, and P. ZoUer, Phys. 
Rev. Lett. 85, 4643 (2000). 

[17] We note that recently there has been established a dy- 
namical connection between perfect fluid cosmologies in 
Einsteinian gravity and the dynamics of expanding Bose- 
Einstein condensates, not related to the kinematical anal- 
ogy considered here; J. E. Lidsey, Class. Quantum Grav. 
21, 777 (2004). 

[18] U. Leonhardt, T. Kiss, and P. Ohberg, J. Opt. B: Quan- 
tum Semiclassical Opt. 5, S42 (2003). 

[19] J.R. Anglin and W. Ketterle, Nature 416, 211 (2002). 

[20] A. E. Leanhardt et al, Science 301, 1513 (2003). 

[21] W. de Sitter, Mon. Not. R. Astron. Soc. 78, 3 (1917). 

[22] A. Recati et al, arXiv:cond-mat/0212413 

[23] P.O. Fedichev and U. R. Fischer, Phys. Rev. Lett. 91, 
240407 (2003). 

[24] Yu. Kagan, E. L. Surkov, and G. V. Shlyapnikov, Phys. 

Rev. A 54, R1753 (1996); Y. Castin and R. Dum, Phys. 

Rev. Lett. 77, 5315 (1996). 
[25] L. P. Pitaevskii and A. Rosch, Phys. Rev. A 55, R853 

(1997). 

[26] High quality factors for breathing mode (monopole ra- 
dial) oscillations in the strongly elongated limit a;||/tJx 
have been obtained experimentally by F. Chevy, V. 
Bretin, P. Rosenbusch, K. W. Madison, and J. Dalibard, 
Phys. Rev. Lett. 88, 250402 (2002). 

[27] G. Baym and C. J. Pethick, Phys. Rev. Lett. 76, 6 (1996). 

[28] E. Zaremba, Phys. Rev. A 57, 518 (1998). 

[29] VI. S. Dotsenko and V. A. Fateev, Nucl. Phys. B 240, 312 
(1984). 

[30] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 

[31] P.O. Fedichev and U. R. Fischer, Phys. Rev. A 69, 



9 



033602 (2004). [33] A. Gorlitz et al, Phys. Rev. Lett. 87, 130402 (2001). 

[32] L. I. Men'shikov and A. N. Pinzul, Usp. Fiz. Nauk 165, 
1077 (1995) [Phys. Usp. 38, 1031 (1995)]. 



